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ABSTRACT 

Dynamic analysis of frames involves very rigorous method which always requires software programming to take 
care of cumbersome calculations that are always involved. An easy manual method has been developed here to take 
care of these rigorous calculations that always require software programming. This method was used to analyze a 
three storey plane rectangular rigid frame. The condensed matrix gotten from this method was compared with that of 
Chopra and the result shows high similarity. Eigen value, eigenvector, shear force and base overturning moment of 
the three storey frame were also determined in the analysis. The results showed that plane rectangular rigid frames 
can easily be analyzed with this method which is less cumbersome and easy to apply. However this method cannot 
easily handle more than ten storey frames because of large matrix calculations involved. 

Keywords: Dynamics, eigenvalues, overturning, eigenvectors and frames. 

1.0 INTRODUCTION 

Galileo was one of the first to deal thoroughly with the concept of acceleration and he founded 
dynamics as a branch of natural philosophy. The close interplay of theory and experiment, 
characteristic of this subject, was founded by Italian scientists. Galileo said mathematics is the 
means to decipher the book of nature. Mathematics seeks to discern the outlines of all possible 
abstract structure. This pure mathematics may be applied to every sort of concrete problem 
(Corradi, 2013).Structural dynamics deals with the structures that are subjected to dynamic 
loading. The loads include windstorm, earthquake loads, people, hurricanes, waves, vibration of 
the ground due to a blast nearby, flooding, tornadoes, vibration effect from compound disc 
player, vibrations effect from moving vehicles, vibration effect from heavy construction 
equipment etc. Any structure can be subjected to dynamic loading. Dynamic analysis can be used 
to find dynamic displacements, Eigen value, eigenvector, shear force and bending moment. 
(Ezeh et al, 2014). The problem of dynamic analysis is the rigorous and cumbersome 
calculations that one has to solve to obtain the Eigen values and eigenvectors. Recently, authors 
have developed software that can handle this cumbersome calculation but there is still need to 
develop manual method which could facilitate learning and manual analyses of these dynamic 
loads. This fall into the main objectives of this study which is the algorithm for the dynamic 
analysis of plane rectangular rigid frame subjected to ground motion. This objectives can be 
achieved through the following methods: 

♦ Matrix stiffness method to analyse plane rectangular rigid frame. 

♦ Method of static condensation which reduces the number of degrees of freedom of the 
plane rectangular rigid frame. 

♦ Household characteristic equation, Newton-Raphson’s method and Gaussian elimination 
method in the vibration analysis of plane rectangular rigid frame subjected to ground 
motion. 

♦ Manually determination of the natural frequency of different mode shape (eigenvalues 
and eigenvectors of the frame subjected to ground motion). 

♦ Solving for dynamic base shear force and dynamic overturning moment. 
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2.0 Materials and Methods 

2.1 Elastic Frame Assumptions 

i. Members of the frame are axially rigid (stiff). Hence, there are only two possible 
displacements per joint or node. 

ii. The element of the frame are concieved as finite (lump) mass system. 

iii. The nodal displacement of the frame are numbered as shown in Figure 1 

iv. The U’s stands for the lateral and rotational displacement 

v. Axial deformation of beams, columns and the effect of axial forces on the stiffness of the 
columns were neglected. This is because the inertial effect associated with rotations and 
vertical displacements of the joints are usually not significant. (Clough and Penzien 1995) 

2.1.1 Numbering of the Frame 

The nodes, elements and the deformation or forces acting on the frame were numbered first from 
the base to the top, having in mind that a frame with N number of lateral deformation will also 
have 2N number of rotational deformation i.e figure 1 has 6 number of lateral deformations and 
6 number of rotational deformations. This leaves every node with three DOFs which are: 

i. Axial deformation 

ii. Lateral deformation 

iii. Rotational deformation. 

With axial deformation neglected, the total deformations at each node reduces to two (lateral 
and rotational deformation). 



Figure 1: Three storey frame showing the mass, lateral and rotational displacement. 


2.2 Method of static condensation elements and global stiffness matrix. 

Roy and Chakrabarty, (2009) derived the element stiffness matrix of an axially rigid line 
elements with two nodes as: 
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( 1 ) 

Equation (1) was used to determine the elements stiffness matrices for every element of the 
frame. After which the various elements stiffness matrices were combined to form the global 
stiffness matrix. 

Where P t and Pj are the initial and final loads applied to the element and Mjand Mj is the initial 
and final moments generated from applied load on the element. V t , Vj and 6 t , Gj are the initial 
and final lateral and rotational deformations generated from the effect of the load respectively. 

Static Condensation. 

This method eliminate those degrees of freedom without masses assigned to them. This happens 
by eliminating the rotational deformations and leaving the structure with only lateral 
deformations. Thereby reducing the global matrix to one-third of its original size (Anderson and 
Naeim, 2012).Thus the global matrix was partitioned as shown in matrix A and the necessary 
condensation were carried out as shown. 


Matrix A = 


Ktt 

K t0 

Kot 

Koo 


k 0 Q and /cj t were determined from K 00 and K ot of matrix A and substituted in equation (2) to 
obtain the condensed matrix. 

^tt = ktt " ^ot^oo^ot (2) 

The steps are summarized below 

i. Partition the global stiffness matrix like the static condesation pattern of matrix A. 

ii. Obtain /cj t = k t0 and k^ 

iii. Obtain T = -koo^ot hence, -/Cot^oo^^ot 

iv. Substitute for k tt = k tt - ko t kook ot 

v. k tt is the condensed stiffness matrix which ensures that the matrix has been reduced to 
one-third of its original size. 

3.0 Vibration analysis using household characteristic equation and Newton- 
Raphsoniteretion method. 

Eigenvalues are obtained by solving for w 2 in equation (3). Thus 

Det(k tt - w 2 m tt ) (3) 
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Where w 2 = A = eigenvalue, m tt = lump mass matrix. 

Substituting the values of condensed matrix fc tt and lump mass matrix m tt in equation (3) and 
solving the determinant of equation (3) by applying Szilard 2004, household characteristic 
equation gives a polynimial equation of the structure. This polynomial equation is solved by 
applying Bird 2010, Newton - Raphsons iterations method for solving root of polynomial 
equation which will yield various eigenvalues of the frame. 

Below are the illustration on how to generate the polynomial. 


a ll~ & 0-12 a 13 . a ln 

a 21 a 22- & a 23 . a 2n 

a 31 a 32 a 33‘ ^. a 3 n 

A = det . (4) 


a nl a n2 a n3 . O-nn ^ 

The household characteristic equation were derived by Szilard, (2004) as expressed. 


/L n + M n_1 + k 2 /L n ~ 2 + ..,+ k n ^A + k n = 0 (5) 

Where coefficients k x , k 2 , . k n are defined by 

k x = -T x (6a) 

k 2 = -\{k x T t + T 2 ) (6b) 

k 3 = -^k 2 T i + k x T 2 + T 3 ) (6c) 

K = ~~ik n - 1 T 1 + k n _ 2 T 2 +.+ kiT^ + T n ) (6d) 

In equation (6), T t , T 2 , T 3 , .T n represent the traces of the matrices A n , which are defined 


as the sum of the diagonal elements of the corresponding matrix. Thus we can write for matrix A 
7\ = trace A 
T 2 = trace A 2 


T n = trace A n 

The trace of a matrix is the summation of the diagonal element in the matrix. Consider matrix A 
of equation (4). 

7\ = trace A = a xl + a 22 + a 33 +.+ a nn (7a) 

T 2 = trace A 2 = a 2 ± + a 22 + a 33 +.+ a£ n (7b) 

T 3 = trace A 3 = + a| 2 + o. 33 +.+ (7c) 


T n = trace A n = a 31 + a 22 + a 33 +.+ a£ n (7d) 
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Substituting the value of equation (7) into equation (6) yielded the value of /c l5 k 2 , k 3 , . k n . 

Thus the value of k n and n for A n are substituted in equation (5) to obtain the polynomial 
equation of n degrees. The polynomial equation is solved by applying Newton - Raphson’s 
iteration as explained. 


3.1 Newton - Raphson’s Iteration 

In Bird (2010), the Newton - Raphson iteration popularly known as the Newton-Raphson’s 
method, is stated as follows; “if ri is the approximate value of a real root of the equation 
f(x) = 0, then a closer approximation r 2 , to the root is given by; 

fQt) 

/ 7 (a) 


r 2 =r i~Z 


(8) 


Equation (8) above is used to solve the polynomial equation of any degree by trial and error. 
Where 

r 2 = root of the polynomial equation gotten by solving equation (8) 
r x = assumed root for trials 


f(r x ) = the value obtained when r x is substituted in the polynomial equation. 

/ / ( r i) = the value obtained when r x is substituted in the differential of the polynomial equation. 
A value for the root is assumed as r x and this value is substituted in the polynomial equation to 
obtain f(?i). The polynomial equation is differentiated and the same r x is substituted in the 
polynomial as the root to obtain / / (r 1 ). The result is then substituted in equation (8) to obtain r 2 
and r 2 is subtstituted back in the polynomial to check if it is a true root of the polynomial 
equation. The true root of a polynomial equation will give zero when substituted in the equation. 
Continue with the same proceedure till all the roots of the polynomial are obtained. 


3.2 Steps to follow when solving for the root of the polynomial equation using equation 

( 8 ). 

i. Choose a trial root value r lf r x must be a positive integer since the eigenvalue can not be 
a negative integer. 

ii. Obtain the differential of the polynomial equation of which its roots are required. 

iii. Substitute the trial root value r x in the polynomial to obtain a value f(r x ) 

iv. Substitute the trial root value r x in the differential of the polynomial to obtain another 
value f’(r : 1 ) 

v. Substitute the value of r 1; f(r x ) and / 7 (r 1 ) in equation (8) to obtain r 2 . 

vi. Check if r 2 will satisfy the polynomial equation by substituting r 2 in the polynomial 
equation to see if the polynomial will give zero. 

vii. If r 2 did not satisfy the polynomial equation, then increase r ± and repeat step 3 to 6. 

viii. Continue step 7 until a value of r 2 that satisfy the polynomial is obtained and this 
value is one of the root of the polynomial depending on the degree of the polynomial. 

ix. Increase r x again and repeat step iii to viii to obtain the second root of the polynomial 
equation. 

x. Continue step 9 until all the roots of the polynomial are obatained (depending on the 
degree of the polynomial). These roots are the required eigenvalues. 
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3.3 How to solve for Eigenvector or Natural Mode <P Using Gaussian Elimination 
Method. 

The eigenvectors are obtained by solving equation (9) below. 

(k t t - 6nm tt )<I> nj = 0 (9) 

The number of eigenvalues in a system determ in e the number of eigenvectors as well. The 
number of eigenvectors in a system are the square of the number of eigenvalues provided the 
number is not one. Each eigenvalues generate eigenvector that equals the number of eigenvalues 
of the system, i.e. For the case of three story frame of figure (1), n = 3 and A 1 will generate 0 ll5 
0 12 &(P 13 , X 2 will generate d> 21 , < & 22& < &23 and X 3 will generate 0 31 , 0 32 &0 33 . From the last 
statement we noted that the system has three eigenvalues and nine eigenvectors. 

Substituting the value of k tt A n &m tt in equation (9) gives a matrix of n degrees. To obtain the 
eigenvectors, Gaussian elimination method is applied to the linear simultaneous equations i.e. 
one of the eigenvectors from the singular matrix has to be assumed to be one and then the 
fraction of others are solved for. 

The set of the eigenvector simultaneus equation can be given as follows. 


an0n + a 12 $ 12 + a 13 0 13 +.+ a ln 0 ln - 0 

a 21^*21 + 0-22^22 + d-23^23 + . T d. 2n 02n ~ 0 

^31*3^31 + a 32 0 32 + CL 33 0 33 T . H~ d 3n 0 3n = 0 


drn&nl + a n2 d> n2 + 0-n3^n3 +.+ O-nn^nn = 0 


-a 11 a 12 a 13 . 

. a ln' 

\ 0U1 ] 


-0- 

a 21 a 22 a 23 . 

. a 2n 

<t>n2 


0 

a 31 a 32 a 33 . 

. a 3n 

^n3 

= 

0 

-CL n i a n2 a n3 . 

. a nn . 

-®nn- 


0 


3.4 Steps to follow when solving for Eigenvectors 

i. Let 0 nn = 1 

ii. Divide all the row by the coefficient of 0 nl 

iii. Subtract every other row from the first row to make the coefficient of 0 nl of every other 
row apart from first row to be zero. 

iv. Divide every other row apart from the first row by the coefficient of 0 n2 to make the 
coefficient of 0 n2 to be one (1). 

v. Subtract every other rows apart from first row from the second row to make their 
coefficient 0 n2 to be zero. 

vi. Divide every other row apart from first and second row by the coefficient of $ n3 to make 
their coefficient 0 n3 to be one (1). 

vii. Subtract every other rows apart from first and second row from the third row to make 
their coefficient <P n3 to be zero. 

viii. Continue in this manner till the last row become solvable. 

ix. Repeat step 1 to 8 for all the eigenvalues (A n , where n = 1, 2, 3, ...) 
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Table 1: Modal Static and Dynamic Response, r 


Modal Static Response, r 

Modal Dynamic Response, r 

Vi = V£=ZUs,n 

Mi == M^ = ^U(hi-hi)s in 

K = On( 0 

Mi = M%A n (t) 

V b == v& = ZU s ln = r„ii5 s M- 
M b == M b n = Zj=ihjSjn = r n Lb = h n M n 

Uj == Ujn = (UOtjn 

*1 = A& = zu (r n /w„ 2 )(0 /n - 4>,. lin ) 

V b = On(t) 

M b = M&A n ( t) 

Uj = uf n A n { t) 

Aj = A%A n ( t) 




Where r„i5 = ^K = ^l‘n = ZJU h,mj 0 jn 


Where V bn and M bn are the base shear and overturning moment due to nth mode. 

3.5 Overall steps to the present study 

STEP1: Number the frame. When numbering the frame, number the joint separately, also 
number the two component members (beams and column) together and then number the 
displacements separately. For each joint, we numbered the lateral and rotational displacement. 

STEP2: Determine the stiffness matrix. The stiffness matrix equation of equation (1) was used to 
determine the stiffness matrix of the frame. Each element has two nodes and every node has two 
deformations, lateral and rotational deformation. But in a beam element, the lateral deformations 
are neglected because it has been taken care of in the neighboring joint column. 

STEP3: Add the stiffness matrices for all the elements to obtain the global stiffness matrix. 

STEP4: Determine the condensed stiffness matrix. The condensed stiffness matrix is determined 
by the method of static condensation. This method eliminates from dynamic analysis those 
degrees of freedom of a structure to which zero mass is assigned. 

STEP5: Determine the lump mass matrix. 

STEP6: Determine the eigenvalue from the frequency equation (3). Solve equation (3) by 
applying Szilard household characteristic equation and Birds (2010) Newton-Ralphson 
procedures for solving root of polynomial equations. 

STEP7: Determine the natural modes or Eigenvectors O n by Gaussian elimination. 

STEP8: Draw the graphical representation of the Eigenvectors. 

STEP9: Determine the base shear forces. 

STEP10: Determine the base overturning moments. 

Chopra 2007, summarised the formular for calculating modal static and modal dynamic 
responses with tables as shown in Table 1. 
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4.0 Conclusions. 

Based on the results of this research, the following conclusions were drawn; 

i. A three storey Plane rectangular rigid frames subjected to ground motions was analyzed. 

ii. The steps marshalled out in this work yielded good results when applied to dynamic 
frame. 

iii. The method developed showed high level of accuracy when the result of the condensed 
matrix obtained with it was compared with that of Chopra. 

iv. The method is an efficient and satisfactory method for dynamic analysis of frame 
subjected to ground motion. 

5.0 Recommendations. 

The following reco mm endations were made: 

i. The develop highly recommends this method in the dynamic analysis of frames subjected 
to ground motion 

ii. Future research work should consider using matrix stiffness method to solve for 
rectangular frames that has consistence mass and compare its convergence with plane 
rectangular frames that has lump mass. 


APPENDICE 1 
Example: 

i. The figures “a” below showed a three storey frame with lumped masses subjected to 
lateral forces, together with its properties; in addition, the flexural rigidities are given 
in El’s for colu mns and beams. Determine the condensed matrix. Neglect the axial 
deformation of the members. Assume the members to be massless and determine the 
eigenvectors and overturning moments. 



Appendix 1: Three storey rectangular plane rigid frame 
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Solution 

Equation (1) was used to determine the elements stiffness matrices for every element of the 
frame. After which the various elements stiffness matrices were combined to form the global 
stiffness matrix thereby forming a nine by nine global matrix as shown. 



Condensed Matrix 

The global matrix was partitioned as shown in matrix A and the necessary condensation were 
carried out using equations (2). The results were shown below. 


40.84645161 

-23.25677 

5.109677 

-23.25677419 

31.091613 

-14.2452 

5.109677419 

-14.24516 

10.06452 


Solving For Eigenvalues. 

The mass matrix and condensed matrix were substituted in equation (3) to obtain the eigenvalues 
as shown below 


Mass Matrix- 


1 

0 

0 

0 

1 

0 

0 

0 

0.5 


Eigenvalue 1 


38.72185161 

-23.25677 

5.109677 

-23.25677419 

28.967013 

-14.2452 

5.109677419 

-14.24516 

9.002216 
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Eigenvalue 2 


28.64015161 

-23.25677 

5.109677 

-23.25677419 

18.885313 

-14.2452 

5.109677419 

-14.24516 

3.961366 


Eigenvalue 3 


18.11095161 

-23.25677 

5.109677 

-23.25677419 

8.3561129 

-14.2452 

5.109677419 

-14.24516 

-1.30323 


Solving for Eigenvectors 

The eigenvectors shown are obtained by solving equation (9). 

Eigenvector Graph _ 

Graph (Y) _ Mode 1 (X) _ Mode 2 (X) _ Mode 3 (X) 

0 0 0 0 

4 0.2708 -0.0142 -1.0115 

7 0.6589 0.2028 -0.5245 

10 1 1 1 


Solving for base shear force and overturning moment 

The modal properties, base shear force due to nth mode and overturning moment due to nth 
mode are obtained with the help of table (1). 


The Modal Properties 


S/S 

Mn 

Ln 


Rn 

Yb 

1 

1.01 

1.43 


1.419 

2.029 

2 

0.54 

0.689 


1.272 

0.876 

3 

1.8 

-1.04 


-0.576 

0.597 

Base shear due to nth mode 

S/S 



Ybbn 



1 



2.03 



2 



0.88 



3 



0.6 




Overturning moment due to nth mode 


S/S 

L 

Mbn 

1 

6.45 

9.15 

2 

3.77 

4.794 

3 

-2.13 

1.229 


Efficacy of the method for condensed matrices 

Chopra (2007), showed the condensed matrix of a three storey plane rectangular frame with 
lumped masses subjected to lateral forces, together with its properties; in addition, the flexural 
rigidity is El for columns and beams. The problem was given in appendix 1 (Example). 


Chopra condensed matrix (matrix C) are listed below 
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1 40.85 -23.26 5.11 ' 

-23.26 31.09 -14.25 

5.11 -14.25 10.06 . 

Present Method condensed matrix (matrix D) are listed below 


' 40.8465 -23.2568 5.10968 ' 

-23.25677 31.0961 -14.2452 

. 5.109677 -14.2452 10.0645 . 


From the result of the condensed matrix obtained, it was observed that Chopra’s solution and the 
present study values for the condensed matrix yielded approximately the same values. 

Thus, the method is satisfactory and can be relied upon in dynamic analysis. 
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